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We analyze the thermodynamical state of nuclear matter in transport 

' calculations of heavy-ion reactions. In particular we determine temperatures 

o ■ 
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and radial flow parameters from an analysis of fragment energy spectra and 

■ compare to local microscopic temperatures obtained from an analysis of local 
momentum space distributions. The analysis shows that the spectator reaches 
an equilibrated freeze-out configuration which undergoes simultaneous frag- 

(<~^ ■ mentation. The fragments from the participant region, on the other hand, do 

■ not seem to come from a common fragmenting source in thermodaynamical 

in 
o 

0^ 
0^ 



X 



equilibrium. 
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2 ■ One of the major interests in the study of intermediate energy heavy-ion collisions is the 



understanding of the multifragmentation phenomenon and its connection with liquid-gas 
phase transitions For this it has to be assumed that in a heavy ion collision at some 
stage a part of the system is both in thermodynamical equilibrium and instable. Such a con- 
figuration is often termed a freeze-out configuration. The multifragmentation process would 
reflect the parameters of this source, i.e. its temperature, density, and perhaps collective 
(radial) flow pattern. Experimentally the fragment energy spectra are described in terms 
of such freeze-out models to extract these parameters. On the other hand, information on 
these quantities is also sought from other observables, in particular from isotope ratios of 
fragments, and from excited state populations. In the past the conclusions drawn from these 
different sources have often been in conflict with each other. 

One way to study whether this szenario is applicable is to analyze the results of transport 
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calculations of heavy ion collisions g-yj. Since such calculations reproduce reasonably well 
the asymptotic observables it should be meanigful to look also into their intermediate time 
behaviour to see whether, where, and when such freeze-out configurations exist. In ref. ^ 
a method was developed to determine thermodynamical variables locally by analyzing the 
local momentum distribution even in the presence of possible anisotropics. 

In the present work we apply this analysis to intermediate energy collisions of Au + Au, 
which were also studied extensively experimentally [^|T0(]. We want to establish whether the 
concept of a freeze-out configuration used in the experimental analysis of fragment kinetic 
energy spectra is supported in realistic transport calculations. Fragments are described 
using a coalescence algorithm as detailed later. We study, in particular, the participant 
and spectator regions which represent rather clean thermodynamical situations in a heavy 
ion collision. In this way we also try to clarify some of the discrepancies between different 
methods of temperature definition. 

We base our investigation on relativistic transport calculations of the Boltzmann- 
Nordheim-Vlasov type. In this work we use, in particular, the relativistic Landau- Vlasov 
approach which was described in detail in ref. It uses Gaussian test particles in co- 
ordinate and momentum space and thus allows to construct locally a smooth momentum 
distribution. For the self energies in the transport calculation we have adopted the non- 



linear parametrization NL2 |jTT|. In ref. |jT2| we compared this parametrization to more 



realistic non-equilibrium self energies based on Dirac-Brueckner calculations. With respect 
to the thermodynamical variables discussed here, we found no essential differences between 
the two models, and thus we use the simpler NL2 here. A similar analysis with respect 
to fragment kinetic energy spectra in central collisions has been performed previously in 
ref. in the framework of non-relativistic transport theory with special emphasis on the 
dependence on assumptions of the equation of state and in-medium cross sections. Here also 
a weak dependence of radial fiow observables on particular choices of the microscopic input 
was found. 

The microscopic determination of a local temperature from the phase space distribution 



was discussed in detail in refs. [0. Thus we only briefly review the procedure here. The 
local momentum distribution obtained from a transport calculation is subjected to a fit in 
terms of covariant hot Fermi-Dirac distributions of the form 

n{x,k,T) = i J (1) 

1 + exp — (/i* — /c*'u^)/T 

with the temperature T, the effective chemical potential /i*(T) and = E* = \J k*^ + m*^. 
For vanishing temperature eq. (|lD includes the limit of a sharp Fermi ellipsoid with /i* = Ep = 



kp + m*2. The local streaming four- velocity is determined from the local 4-current j'^ 



as = j^/po, where po = yJf^j'^ is the local invariant density. Then the temperature T 
is the only fit parameter to be directly determined from the phase space distribution. In 
this procedure the effect of the potential energy is taken into account by way of effective 
masses and momenta and the Fermi motion of the correct density by the chemical potential 
fi*. Thus this temperature is a local thermodynamic temperature, which in the following we 
denote as Tioc- 

Expression is appropriate for a system in local equilibrium. In a heavy ion collision 
this is generally not the case and would lead to an interpretation of collective kinetic energies 
in terms of temperature. To account for anisotropy effects, e.g. in ref. [|T3| longitudinal and 



perpendicular temperatures have been introduced. In our approach we model anisotropic 
momentum distributions by counter-streaming or 'colliding' nuclear matter i.e., 
by a superposition of two Fermi distributions n^^'^^ = n^^^ + n'-^-' — Sn^^'^\ where Sn^^"^^ = 
V n(^) ■ ?7,(2) guarantees the validity of the Pauli principle and provides a smooth transition 
to one equilibrated system. In it has been demonstrated that this ansatz allows a 

reliable description of the participant and spectator matter at each stage of the reaction. 

Experimentally much of the information about the thermodynamical behaviour in heavy 
ion collisions originates from the analysis of fragment observables. Thus also in the present 
analysis we will need to generate and analyze fragments. The correct and practical procedure 
how to properly decribe fragment production is still very much debated [|T^,0. Here we 
do not enter into this debate but use the simplest algorithm, namely a coalescence model. 



as we have done and decribed it in ref. [0. In brief, we apply phase space coalescence, 
i.e. nucleons form a fragment, if their positions and momenta {xi,pi) satisfy \xi — Xf \ < 
and \pi — Pf\ < Pc- Rc,Pc are parameters which are fitted to reproduce the observed mass 
distributions and thus guarantee a good overall description of the fragment multiplicities. 

Fragment kinetic energy spectra have been analyzed experimentally in the Siemens- 
Rassmussen or blast model [|,^|l8l. In this model the kinetic energies are interpreted in 



terms of a thermalized freeze-out configuration, characterized by a common temperature 
and a radial fiow, i.e. by an isotropically expanding source. The kinetic energies are given 
by 



— ~ y /3^d/3ni/3) eMlE/T) x [-— (^7 + - j - _ cosh a] , (2) 

where p and E are the center of mass momentum and the total energy of the particle with 
mass m, respectively, and where 7"^ = 1-/3^ and a = 'jPp/T. Various assumptions 
have been made for the fiow profile n{i3). A good parametrization is a Fermi-type function 
| |1U| |. However, the results are not very different when using a single fiow velocity, i.e. 
n{f3) ~ 6{(3 — f3f), which we also use here for simplicity. One then has two parameters in 
the fit, namely f3f and the temperature parameter in eq.(^, which we call Tgiope- It is, of 
course, not obvious that Tgiope represents a thermodynamical temperature. One of the aims 
of this investigation is, in fact, to find its significance. The expression (0) has been applied to 
kinetic energy spectra of all fragment masses simultaneously, yielding a global Tgiopeiglobal) , 
or to each fragment mass separately, giving Tsiopei^f). If a global description was achieved, 
it was concluded that a freeze-out configuration exists. We will also test this procedure. 

In the following we apply the above methods to central (6 = fm) and semi-central 
{b = 4.5 fm) Au + Au collisions at Eheam = 0.25 — 0.8 AGeV. This reactions have been 
studied extensively by the ALADIN [p], p!9| , p0| and EOS collaborations with respect to 
temperature and phase transitions. In ref. [0 we have previously studied this reaction at 
one energy, 600 MeV/A, only with respect to local temperatures and thermodynamical 
instabilities. In this work we perform the fragment analysis with the blast model to extract 



and compare slope temperatures and we discuss a wider range of incident energies [^^ . 

The spectator is that part of the system which has not colhded with the other nucleus, 
but which is nevertheless excited due to the shearing-off of part of the nucleus and due 
to absorption of participant particles. In the calculation it is identified as those particles 
which have approximately beam rapidity. It was seen in ref. [0 that it represents a well 
equilibrated piece of nuclear matter at finite temperature. 

In fig. 1 we show the evolution with time of the local temperature and the density 
for the spectator at various incident energies. After the time when the spectator is fully 
developed the properties are rather independent of incident energy which supports the freeze- 
out picture. Also after this time the density and temperature remain rather constant for 
several tens of fm/c, making it an ideal system in order to study the thermodynamical 
evolution of low-density, finite temperature nuclear matter. In ref. we also determined 
pressure and studied the dependence of pressure on density. We found that after about 45 
fm/c the effective compressibility K ~ dP/dplx became negative indicating that the system 
enters a region of spinodal instability and should subsequently break up into fragments. 
At this time we find densities of about p ~ 0.4 — 0.5po and T ~ 5 — 6 MeV, which is 
in good agreement with findings of the ALADIN group based on isotope thermometers 
19| , |20|| . Recently the ALADIN group has also determined kinetic energy spectra of spectator 



fragments and has extracted slope temperatures using eq.(^. It was found that these are 
typically 10 to 12 MeV higher than those measured with the isotope thermometer. 

Applying the coalescence model to the spectator we obtain kinetic energy spectra as 
shown in fig. 2 at 600A MeV for nucleons {Aj = 1) and for fragments with Af > 2 separately. 
We fitted these spectra with the model of eq. (0) in the rest frame of the spectator = 0). 
The Af > 2 spectrum is well fitted with a temperature of T^/ope ~ (17 ± 2) MeV. The 
nucleon spectrum, on the other hand, shows a two-component structure, as was also observed 



experimentally in ref. It is dominated by a low energy component with Tsiope,iow = 

(7.3 ± 3.5) MeV. The high energy component has rather poor statistics in our calculation 
and we interpret it as nucleons from the participant that have entered the spectator region. 



The slope temperature of the low energy component is close to the local temperature Tioc = 
(5 — 6) MeV as discussed above with respect to fig. 1. Thus for nucleons both methods 
of temperature determination consistently are seen to yield about the same result. In fact, 
they should not neccessarily be identical, since Tioc is determined fitting the momentum 
distribution by a Fermi function while Tgiope in eq. is based on a Maxwell-Boltzmann 
distribution, which are not the same at such low temperatures. 

On the other hand the slope temperatures of the fragments are considerably higher than 
those of the nucleons. In fig. 3 we show the slope temperatures separately for the different 
fragment masses and also the local temperature for comparison. There is a rapid increase 
of Tsiope with fragment mass which saturates for Af > 3 around Tgiope ~ 17 MeV, which 



was the temperature determined in fig. 2. The experimental values from ALADIN pO[ also 
shown in fig. 3. were obtained by analogous blast model fits to the measured spectra. It 
can be seen that the slope temperatures from the theoretical calculations and from the data 
agree extremely well. Also the corresponding kinetic energies which range from from 23.7 
MeV (Af = 2) to 28.1 MeV (Af = 8) are in good agreement with the ALADIN data. 

At first sight it is surprising that Tgiope for nucleons and fragments differ from each other 
and also from Tioc- The difference has been interpreted in ref. in terms of the Goldhaber 
model [^, as it has been applied to fragmentation by Bauer [^. When a system of 



fermions of given density and temperature suddenly breaks up the fragment momenta are 
approximately given by the sum of the momenta of the nucleons before the decay. For 
heavier fragments the addition of momenta can be considered as a stochastic process which 
via the central limit theorem leads to Gaussian energy distributions which resemble Maxwell 
distributions and thus contribute to the slope temperature. As discussed by Bauer and also 



in ref. I^Oj this effect increases the slope temperatures by an amount which is of the order 
of the difference between the isotope and the slope temperatures. 

We wanted to see whether a similar effect can explain the mass dependence seen in 
fig. 3. We therefore initialized statistically a system of the mass and temperature of the 
spectator, and subjected it to the same fragmentation procedure (coalescence) and to the 



same fit by eq. (0) as we did for tlie lieavy ion collision. These slope temperatures obtained 
from the statistical model are given in fig. 3 as a band, which corresponds to initializations 
between p = O.Spo and T = 6 MeV and p = 0.4po and T = 5.5 MeV, which cover the 
range of values in fig. 1. It is seen that the model qualitatively explains the increase in 
the slope temperature relative to the local temperature and the increase with fragment 



mass relative to that for nucleons. A similar conclusion was drawn in ref pO[ using the 



results from ref. |2^. This shows that Tgiope is not a thermodynamic temperature. The 



difference relative to the thermodynamic temperature can be understood from the fact, that 
to form a fragment the internal kinetic energies of the nucleons are limited by the coalescence 
condition. Since on the average all nucleons have the same momenta, this means that the 
collective momentum per nucleon of the fragment increases relative to the average. This 
simulates a higher temperature. This effect has been called "contribution of Fermi motion 
to the temperature" . The purpose of using the Goldhaber model here was to demonstrate 
this effect. Whether a Goldhaber model applies to heavy ion collisions, can, of course, be 
debated, but our results are independent of this question. 

Thus we seem to understand fairly well the kinetic energy spectra of the spectator frag- 
ments and we now turn to the participant region. The participant zone in a heavy ion 
collision constitutes another limiting, but still simple case for the investigation of the ther- 
modynamical behaviour of nuclear matter. In contrast to the spectator zone one expects 
a compression- decompression cycle and thus richer phenomena with respect to fragmenta- 
tion. The situation becomes particularly simple if we look at central collisions of symmetric 
systems which experimentally are selected using transverse energy distributions, charged 



particle multiplicities or polar angles near mid-rapidity [p|, p!0| , p4 . 

We begin to characterize the calculated evolution of a collision for the case of Au + Au 
at 600A MeV. A very well developed radial flow pattern appears after about 20 fm/c in 
agreement with findings of other groups 0^,^]. The pressure in this reaction becomes 
isotropic at about 35 fm/c indicating equilibriation. The number of collisions drops to small 
values at about 40 fm/c. This condition we shall call (nucleon) freeze-out. Thus equilibration 



and freeze-out occur rather simultaneously. We find a density at this stage of about normal 
nuclear density and a (local) temperature of about Ti^c ~ 15 MeV in the mid-plane of the 
reaction. 

We now also apply the blast model of eq. (|^) to fragment spectra generated in the 
coalescence model at the end of the collision at about 90 fm/c. The results for the slope 
temperature Tsiope and the mean velocity /?/ are shown in fig. 4 for a common fit to all frag- 
ments with Af > 2 for different incident energies. These are compared to the corresponding 
values extracted by the EOS [|| and FOPI ||9|,p!0[| collaborations by analogous blast model 



fits to charged particle spectra. Our results are in good agreement with the temperatures 
determined by the FOPI collaboration and somewhat lower than those from EOS [§, 
in particular at higher incident energies. For the radial flow, the situation tends to be in 
reverse, in particular with respect to the FOPI results of ref . . This is generally consistent 
with the findings of other groups. E.g. in ref. in a similar approach results were obtained 
for the EOS data, which are close to the data for T^/ope and and above the data for jSf. 
One should keep in mind, however, that Tgiope and /?/ are not independent fit parameters. 
Within the uncertainities of the description by blast model fits there is qualitative agreement 
between calculation and experiment. 

As was done for the spectator we also apply the blast model separately for different 
fragment masses Af. This is shown in fig. 5 at 600 AMeV in the left column. We observe 
that slope temperatures rise and fiow velocities fall with fragment mass in contrast to the 
behaviour for the spectator fragments in fig. 3 where T^iope was about constant. A similar 
behaviour has been seen experimentally at 1 A.GeV in ref. and in calculations in ref. 0. 



It can also be deduced from fragment spectra at 250 A. MeV shown in ref. |T^, which yield 
values very close to the ones given here. 

The fragment mass dependence of Tgiope is thus much stronger and qualitatively different 
than that for the spectator seen in fig. 3. Thus this behaviour cannot be interpreted 
as fragments originating from a common freeze-out configuration, i.e. from a fragmenting 
source. To arrive at an interpretation we have shown on the right column of fig. 5 the local 
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temperatures and flow velocities for different times before the nucleon freeze-out, i.e. for 
t' = t freeze-out — t, with t freeze-out ~ 35/m/c. (We recall that the values at the left of fig. 5 
are obtained at the end of the reaction, i.e. about 90 fm/c.) It is seen that for Af = 1 the 
values at freeze-out are close to the blast model ones, as required. However, for fragment 
masses Af > 1 the slope temperatures and velocities behave qualitatively very similar to 
the local temperatures and fiow velocities at earlier times . This would suggest to interpret 
the fragment temperatures and velocities as signifying that heavier fragments originate at 
times earlier than the nucleon freeze-out. This may not be unreasonable since in order to 
make a heavier fragment one needs higher densities which occur at earlier times and hence 
higher temperatures. However, this does not neccessarily imply that the fragments are really 
formed at this time, since fragments could hardly survive such high temperatures, as also 
discussed in ref. But it could mean that these fragments carry information about this 
stage of the collision. In any case it means that in the participant region fragments are not 
formed in a common equlibrated freeze-out configuration, and that in such a situation slope 
temperatures have to be interpreted with great caution. 

In summary fragmentation phenomena in heavy ion collisions are studied as a means to 
explore the phase diagram of hadronic matter. For this it is neccessary to determine the 
thermodynamical properties of the fragmenting source. One way to do this experimentally 
is to investigate fragment kinetic energy spectra. In theoretical simulations the thermody- 
namical state can be obtained locally in space and time from the phase space distribution. 
In this work we have compared this with the information obtained from the generated frag- 
ment spectra. We apply this method to the spectator and participant regions of relativistic 
Au + Aw-collisions. We find that the spectator represents a well developed, equihbrated 
and instable fragmenting source. The difference in temperature determined from the local 
momentum space (or experimentally from the isotope ratios) and from the kinetic energy 
spectra can be attributed to the Fermi motion in the fragmenting source as discussed in a 
Goldhaber model. In the participant region the local temperature at the nucleon freeze-out 
and the slope temperature from fragment spectra are different from those of the spectator. 
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The slope temperatures rise with fragment mass which might indicate that the fragments 
are not formed in a common, equilibrated source. These investigations should be continued 
using more dynamic methods of fragment formation. 
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FIG. 1. Local temperature (top) and density (bottom) evolution in the spectator in semi-central 
Au + Au reactions (b=4.5 fm) at different beam energies as a function of time. 
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FIG. 2. Energy spectra for nucleons {Af = 1) and fragments {Af > 2) (filled diamonds and 
open circles, respectively) from the spectator for Au + Au at Ei,eam = 600 AMeV. Blast model fits 
are shown (solid curves) which determine slope temperatures. 
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FIG. 3. Spectator slope temperatures for different fragment masses Af for the reaction as in fig. 
^ (diamonds). Also shown is the nucleon local temperature (circle) and the temperature obtained 
from a statistical model (triangle and gray band, see text). 
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FIG. 4. Slope temperatures (top) and radial flow velocities (bottom) from blast model fits to 
fragment {Af > 1) energy spectra for central collisions for different beam energies (filled dots). 
They are compared to data from [^10|. 
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FIG. 5. Slope temperatures (upper row) and flow velocities (lower row) for the same reaction 
as in fig. 4 at Ejj^am = 0.6 AGeV. In the left column for blast model fits for different fragment 
masses and also for A/ > 1 and for all fragments; in the right column the local values from the 
momentum distribution at times before the freeze-out (see text). 
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